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Networks 

Ebrahim Moradi Shahrivar, Mohammad Pirani and Shreyas Sundaram 


Abstract 

We investigate certain structural properties of random interdependent networks. We start by studying a property 
known as r-robustness, which is a strong indicator of the ability of a network to tolerate structural perturbations and 
dynamical attacks. We show that random fc-partite graphs exhibit a threshold for r-robustness, and that this threshold 
is the same as the one for the graph to have minimum degree r. We then extend this characterization to random 
interdependent networks with arbitrary intra-layer topologies. Finally, we characterize the algebraic connectivity of 
such networks, and provide an asymptotically tight rate of growth of this quantity for a certain range of inter¬ 
layer edge formation probabilities. Our results arise from a characterization of the isoperimetric constant of random 
interdependent networks, and yield new insights into the structure and robustness properties of such networks. 

Index Terms 

Random Interdependent Networks, Robustness, Algebraic Connectivity, Isoperimetric Constant. 


I. Introduction 

There is an increasing realization that many large-scale networks are really “networks-of-networks,” consisting 
of interdependencies between different subnetworks (e.g., coupled cyber and physical networks) ||T[-[[6l. Due to the 
prevalence of such networks, their robustness to intentional disruption or natural malfunctions has started to attract 
attention by a variety of researchers |6)-(9). In this paper, we contribute to the understanding of interdependent 
networks by studying the graph-theoretic properties of r-robustness and algebraic connectivity in such networks. 
As we will describe further in the next section, ?’-robustness has strong connotations for the ability of networks to 
withstand structural and dynamical disruptions: it guarantees that the network will remain connected even if up to 
r — 1 nodes are removed from the neighborhood of every node in the network, and facilitates certain consensus 
dynamics that are resilient to adversarial nodes Ifl0l - lfl4l . The algebraic connectivity of a network is the second 
smallest eigenvalue of the Laplacian of that network, and plays a key role in the speed of certain diffusion dynamics 

IE||. 

We focus our analysis on a class of random interdependent networks consisting of k layers (or subnetworks), 
where each edge between nodes in different layers is present independently with certain probability p. Our model 
is fairly general in that we make no assumption on the topologies within the layers, and captures Erdos-Renyi 
graphs and random /c-partite graphs as special cases. We identify a bound p r for the probability of inter-layer 
edge formation p such that for p > p r , random interdependent networks with arbitrary intra-layer topologies are 
guaranteed to be r-robust asymptotically almost surely. For the special case of k- partite random graphs, we show 
that this p r is tight (i.e., it forms a threshold for the property of r-robustness), and furthermore, is also the threshold 
for the minimum degree of the network to be r. This is a potentially surprising result, given that r-robustness is 
a significantly stronger graph property than r-minimum-degree. Recent work has shown that these properties also 
share thresholds in Erdos-Renyi random graphs lfl3l and random intersection graphs ifTTll . and our work in this 
paper adds random k- partite graphs to this list. 

Next, we show that when the inter-layer edge formation probability p satisfies a certain condition, both the 
robustness parameter and the algebraic connectivity of the network grow as ®{np) asymptotically almost surely 
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(where n is the number of nodes in each layer), regardless of the topologies within the layers. Given the key 
role of algebraic connectivity in the speed of consensus dynamics on networks lfT5l . our analysis demonstrates 
the importance of the edges that connect different communities in the network in terms of facilitating information 
spreading, in line with classical findings in the sociology literature lfT6l . Our result on algebraic connectivity 
of random interdependent networks is also applicable to the stochastic block model or planted partition model 
that has been widely studied in the machine learning literature ifTTIl - GOll . While we consider arbitrary intra-layer 
topologies, in the planted partition model it is assumed that the intra-layer edges are also placed randomly with a 
certain probability. Furthermore, the lower bound that we obtain here for A 2 (L) is tighter than the lower bounds 
obtained in the existing planted partition literature for the range of edge formation probabilities that we consider 
m 3 , m3- Both our robustness and algebraic connectivity bounds arise from a characterization that we provide of 
the isoperimetric constant of random interdependent networks. 


II. Graph Definitions and Background 


An undirected graph is denoted by G = {V. E) where V is the set of vertices (or nodes) and E C V x V is the 
set of edges. We denote the set A/j = {vj € V \ ( Vi,Vj ) € E} as the neighbors of node v % £ V' in graph G. The 
degree of node Vi is di = |A/i|, and d rmn and d max are the minimum and maximum degrees of the nodes in the 
graph, respectively. A graph G' = ( V, E') is called a subgraph of G = ( V., E), denoted as G' C G, if V C V and 
E' C E n { V x V'}. For an integer k € Z> 2 , a graph G is A;-partite if its vertex set can be partitioned into k sets 
Vi, V 2 ,.. •, 14 suc h that there are no edges between nodes within any of those sets. 

The edge-boundary of a set of nodes S C V is given by dS = {(v t , Vj) £ E \ V{ £ S,Vj £ h \ S}. The 
isoperimetric constant of G is defined as itTTI 


i(G) = min 

AcV,\A\<% 


m 
w' 


(i) 


By choosing A as the vertex with the smallest degree we obtain i(G) < d nnn . 

The adjacency matrix for the graph is a matrix A £ {0,l} nxn whose (/', j) entry is 1 if ( Vi,Vj ) £ £, and 
zero otherwise. The Laplacian matrix for the graph is given by L = D — A, where D is the degree matrix with 
D = diag(d], do,..., d n ). For an undirected graph G, the Laplacian L is a symmetric matrix with real eigenvalues 
that can be ordered sequentially as 0 = Ai(L) < A 2 (L) < ■ ■ ■ < A n (L) < 2d max - The second smallest eigenvalue 
A 2 (L) is termed the algebraic connectivity of the graph and satisfies the bound ifTTll 


i(Gf 


< A 2 (L) < 2i(G). 


Finally, we will use the following consequence of the Cauchy-Schwartz inequality: 


( 2 ) 


E' 

i=1 


> 


zL A 


( 3 ) 


where s* £ M for 1 < i < k. 


A. r-Robustness of Networks 

Early work on the robustness of networks to structural and dynamical disruptions focused on the notion of 
graph-connectivity, defined as the smallest number of nodes that have to be removed to disconnect the network 
ll22l . A network is said to be r-connected if it has connectivity at least r. In this paper, we will consider a stronger 
metric known as r-robustness, given by the following definition. 

Definition 1 ( lUW): Let r £ N. A subset S of nodes in the graph G = (V,E) is said to be r-reachable if 
there exists a node V{ £ S such that A) \ S\ > r. A graph G = ( V , E) is said to be r-robust if for every pair of 
nonempty, disjoint subsets of V, at least one of them is r-reachable. □ 

Simply put, an r-reachable set contains a node that has r neighbors outside that set, and an r-robust graph has 
the property that no matter how one chooses two disjoint nonempty sets, at least one of those sets is r-reachable. 
This notion carries the following important implications: 

• If network G is r-robust, then it is at least r-connected and has minimum degree at least r fTOl l. 
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• An r-robust network remains connected even after removing up to r — 1 nodes from the neighborhood of every 
remaining node fl3l . 

• Consider the following class of consensus dynamics where each node starts with a scalar real value. At each 
iteration, it discards the highest F and lowest F values in its neighborhood (for some F e N), and updates 
its value as a weighted average of the remaining values. It was shown in ifTOll . lfl2l that if the network is 
(2F + l)-robust, all nodes that follow these dynamics will reach consensus even if there are up to F arbitrarily 
behaving malicious nodes in the neighborhood of every normal node. 

Thus, r-robustness is a stronger property than r-minimum-degree and r-connectivity. Indeed, the gap between 
the robustness and connectivity (and minimum degree) parameters can be arbitrarily large, as illustrated by the 
bipartite graph shown in Fig. |l(a)| That graph has minimum degree n/4 and connectivity n/4. However, if we 
consider the disjoint subsets V\ U V 2 and V 3 U V 4 , neither one of those sets contains a node that has more than 1 
neighbor outside its own set. Thus, this graph is only 1-robust. 

The following result shows that the isoperimetric constant i(G) defined in <[T|) provides a lower bound on the 
robustness parameter. 

Lemma 1: Let r be a positive integer. If i(G) > r — 1, then the graph is at least r-robust. 

Proof: If i(G) > r — 1, then every set of nodes S C V of size up to has at least (r — l)|Sj + 1 edges leaving 
that set (by the definition of i(G)). By the pigeonhole principle, at least one node in S has at least r neighbors 
outside S. Now consider any two disjoint non-empty sets Si and ,S'y . At least one of these sets has size at most (j, 
and thus is r-reachable. Therefore, the graph is r-robust. ■ 

As an example of Lemma [T] consider the graph shown in Fig. |l(a)| Graph G has isoperimetric constant of at 
most 0.5 (since the edge boundary of V 1 UV 2 has size n/4), but is 1-robust. It is worth mentioning that it is possible 
to construct a graph such that difference between the robustness parameter and isoperimetric constant is arbitrarily 
large. For instance, given any arbitrary t € N and n sufficiently large, assume that the interconnection topology 
between partitions Vy and V 3 in graph G is f-regular (i.e., each node in V 2 (V 3 ) is connected to exactly t nodes in 
V 3 (Vy)) and the rest of the graph is the same as the structure shown in Fig. |l(a)| Then the isoperimetric constant 
of the graph G is at most i/2 by considering the set V\ U Hy However, one can show that G is f-robust. 

We have summarized the relationships between these different graph-theoretic measures of robustness in Fig. 1 1 (b)| 
To give additional context to the properties highlighted above, we note the following. Consider an arbitrary partition 
of the nodes of a graph into two sets. An r-connected graph guarantees that the nodes in one of the sets collectively 
have r neighbors outside that set. An r-robust graph guarantees that there is a node in one of the sets that by 
itself has r neighbors outside that set. A graph with i(G) > r guarantees that each node in one of the sets has r 
neighbors outside that set on average. 

In the rest of the paper, we will study random interdependent networks and show how these various properties 
are related in such networks. 


/' /' /' /' 

/ ' < ' / ' , ' 



(a) 


G is a graph on n nodes 

dmin(G') = T 

G is r-connected 


G is r-robust 



i(G ) > r — 1 



(b) 


Fig. 1: (a) Graph G = (V, E) with V = V\ U V 2 U V3 U V4 and |Vj| = 1 < i < 4. All of the nodes in Vj (V3) are 

connected to all of the nodes in V2 (V4). Furthermore, there is a one to one connection between nodes in V 2 and nodes in 
V 3 . (b) Relationships between different notions of robustness. 
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III. Random Interdependent Networks 

We start by formally defining the notion of (random) interdependent networks that we consider in this paper. 

Definition 2: An interdependent network G is denoted by a tuple G = (G i, Go ,..., G . G p ) where G t = ( V t , Efi 
for f = 1,2.... . k are called the layers of the network G, and G p = {V\ U V 2 U ... U14, E p ) is a /.'-partite network 
with E p C U i^jVi x Vj specifying the interconnection (or inter-layer) topology. □ 

For the rest of this paper, we assume that \V\ | = |V 2 1 = • ■ ■ = |I4| = n and the number of layers k is at least 2. 

Definition 3: Define the probability space (Q n ,E n ,F n ), where the sample space Q„ consists of all possible inter¬ 
dependent networks (Gi, G 2 ,..., G}.. G p ) and the index n € N denotes the number of nodes in each layer. The cr- 
algebra T n is the power set of Q n and the probability measure P n associates a probability P(Gi, G 2 ,..., G*, G p ) to 
each network G = (G 1 , G 2 , ■ ■ ■, Gk, G p ). A random interdependent network is a network G = (G\. G 2 , ■ ■ ■, 6 '/., G p ) 
drawn from il n according to the given probability distribution. □ 

Note that deterministic structures for the layers or interconnections can be obtained as a special case of the above 
definition where P(Gi, G 2 , ■ ■ ■ ,Gk, G p ) is 0 for interdependent networks not containing those specific structures; 
for instance, a random k-partite graph is obtained by allocating a probability of 0 to interdependent networks where 
any of the Gi for 1 < i < k is nonempty. Through an abuse of notation, we will refer to random A:-partite graphs 
simply by G p in this paper. Similarly, Erdos-Renyi random graphs on kn nodes are obtained as a special case of 
the above definition by choosing the edges in G\ , G 2 , ■ ■ ■, Gk and G p independently with a common probability p. 

In this paper, we will focus on the case where G p is independent of G t for 1 < i < k. Specifically, we assume 
that each possible edge of the A:-partite network G p is present independently with probability p (which can be a 
function of n). We refer to this as a random interdependent network with Bernoulli interconnections. We will be 
characterizing certain properties of such networks as n gets large, captured by the following definition. 

Definition 4: For a random interdependent network, we say that a property V holds asymptotically almost surely 
(a.a.s.) if the probability measure of the set of graphs with property V (over the probability space (D n , _F n , P n )) 
tends to 1 as n A 00 . □ 


IV. Robustness of Random interdependent networks 


We will first consider random /c-partite networks (a special case of random interdependent networks, as explained 
in the previous section), and show that they exhibit phase transitions at certain thresholds for the probability p, 
defined as follows. 

Definition 5: Consider a function t(n) = with g(n ) -> 00 as n A 00 , and a function x = o{g{n)) which 
satisfies x -A 00 as n —>• 00 . Then t(n) is said to be a (sharp) threshold function for a graph property V if 

1 ) property V a.a.s. holds when pin) = g ( n ^ +x , and 

2 ) property V a.a.s. does not hold for p[n) = 

□ 

Loosely speaking, if the probability of adding an edge between layers is “larger” than the threshold t(n), then 
property V holds; if p is “smaller” than the threshold t(n). then V does not hold. The following theorem is one 
of our main results and characterizes the threshold for r -robustness (and r-minimum-degree) of random /;:-partitc 
networks. 

Theorem 1: For any positive integers r and k > 2, 

In n + (?’ — 1 ) In In n 
t{n) = - 


is a threshold for r-robustness of random fc-partite graphs with Bernoulli interconnections. It is also a threshold for 
the fc-partite network to have a minimum degree r. 

In order to prove this theorem, we show the following stronger result (see Fig. 1(b) 1 that i(G) > r — 1 when the 
probability of edge formation is above the given threshold. 

Lemma 2: Consider a random A:-partite graph G p = (V] U Ip U • • • U 14 ,E p ) with node sets V) = {(i — l)n + 
1, (i — l)n + 2 for 1 < i < k. Assume that the edge formation probability p = p{n) is such that 

p(n) = lun +h^ where r £ N is a constant such that r > 1 and x = x(n) is some function satisfying 

x = o(lnlnn) and x oc as n oc. Define S r as the property that for any set of nodes with size of at most 
|_A;n/2j, we have |ckS'| > (r — 1)|S|. Then G p has property S r a.a.s. In other words, we have i(G p ) > r — 1 a.a.s. 
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The proof of the above lemma is provided in Appendix [A] The idea behind the proof is to first use a union-bound 
to upper bound the probability that for some set S of size [kn/2\ or less, we have |cfS| < | S\(r — 1). Via the use 
of algebraic manipulations and upper bounds, we then show that this probability goes to zero when the inter-edge 
formation probability p(n) satisfies the conditions given in the lemma. 

Using the above lemma, we can prove Theorem Q] 

Proof of Theorem [7} Consider a random A;-partite graph G p of the form described in the theorem with edge 
formation probability p(n) = —" + i)^ lnn+X » where r € N is a constant and x = x(n) is some function satisfying 
x = o(lnlnn) and x —>• oo as n —>• oo. By Lemma O we know that i(G) > r — 1. Therefore, by Lemma [I] the 
random fc-partite graph G p is at least r-robust a.a.s. This also implies that it has minimum degree at least r a.a.s. 
(by the relationships shown in Fig. 1 1 (b)| ). 

Next we have to show that for p(n) = — " + ^' 111 n ~ x , where x = x(n) is some function satisfying x = 
o(lnlnn) and x —>• oo as n —>• oo, a random A-partitc graph is asymptotically almost surely not r-robust. Let 
G p = (Vi U Vo U • • • U 14, E p ) be a random A:-partitc graph. Consider the vertex set V\ = { v\ ,..., v n }, and define 
the random variable A = X\ + X 2 + • • • + X n where A, = I if the degree of node v t is less than r and zero 
otherwise. The goal is to show that if p{n) = ( I ^ lnn ~ x , then Pr(A = 0) —>• 0 asymptotically almost surely. 

This means that for the specified p(n), there exists a node in V] with degree less than r with high probability. 
Since any r-robust graph must have minimum degree of at least r, we will have the required result. 

The random variable X is zero if and only if Aj = 0 for 1 < i < n. The random variables A, and A j are 
identically distributed and independent when i f j and thus we have 

Pr(A = 0) = Pr(Ai = 0) n = (1 - Pr(Ai. = l)) n < e - nPr ( Xl=1) , (4) 


where the last inequality is due to the fact that 1 — p < e p for p > 0. Now, note that 

r— 1 

nPr(Ai = 1) = 

i =0 

'n(k — 1 ) N 


n(k — 1 ) 
i 


> n 


> n 


r — 1 
n(k — 1 ) 
r — 1 


pfl-p^-V-i 

p r ~ 1 (l — pj n ( k ~ 1 )- r + 1 

p r ~ 1 (^~P) n(k ~ 1 \ 

r— 1 


( 5 ) 


where the last inequality is obtained by using the fact that 0 < (1 — p) r ~ l < 1 for r > 1. Using the fact that 
( n{ r fc _ 7 i 1) ) = n (n^ 1 ) for constant r and k, and (1 — = e w(fe-t)in(i-p) _ f]( e ~ ri ( fc ^ 1 )p^ when np 2 —> 0 

(which is satisfied for the function p that we are considering above), the inequality © becomes 


nPr(Ai = 1) = U yn r p‘ 
Substituting p = ln " +( '.J 1 ^ 1 , 11111 " -t and simplifying, we obtain 


j r n r ~ l p- n (k-T)P 


(k—l)n 

nPr(Ai = 1) = Tl 


(In n + (r — 1 ) In In n — x) 


r —1 


e x = n(e x ). 


(In n) r ~ l 

Thus we must have that linin^oo nPr(Ai = 1) = 00 , which proves that Pr(A = 0) —>• 0 as n —> 00 (from ©). 
Therefore, there exists a node with degree less than r and thus the random /c-partite graph G p is not r-robust for 

/ n _ Inrt+(r—1) In Inn—q: 

P\ n > ~ (fc-l)n ' ■ 

As described in the introduction, the above result indicates that the properties of r-robustness and r-minimum- 

degree (and correspondingly, r-connectivity) all share the same threshold function in random /c-partite graphs, 
despite the fact that r-robustness is a significantly stronger property than the other two properties. In particular, 
this indicates that above the given threshold, random A;-partite networks possess stronger robustness properties than 
simply being r-connected: they can withstand the removal of a large number of nodes (up to r — 1 from every 
neighborhood), and facilitate purely local consensus dynamics that are resilient to a large number of malicious 
nodes (up to i n the neighborhood of every normal node). 










6 


A. General Random Interdependent Networks 

With the sharp threshold given by Theorem |T| for random /. -partite graphs in hand, we now consider general 
random interdependent networks with arbitrary topologies within the layers. Note that any general random interde¬ 
pendent network can be obtained by first drawing a random /.--partite graph, and then adding additional edges to 
fill out the layers. Using the fact that r-robustness is a monotonic graph property (i.e., adding edges to an / -robust 
graph does not decrease the robustness parameter), we obtain the following result. 

Corollary 1: Consider a random interdependent graph G = (G\. G 2 , ..., G}., G p ) with Bernoulli interconnections. 
Assume that the inter-layer edge formation probability satisfies p(n) > lnn + n+x , r € Z>i and x = x(n) 
is some function satisfying x = o(lnlnn) and x 00 as n —>• 00 . Then G is r-robust. □ 

The above result shows that if the inter-layer edge formation probability between all layers of the /’-layer graph 
G is higher than the threshold for r-robustness of a k- partite network, then G is a r-robust network, regardless of 
the probability distribution over the topologies within the layers. 


B. Unbounded Robustness in Random Interdependent Networks 

The previous results (Theorem Q] in particular) established inter-layer edge formation probabilities that cause 
random interdependent networks to be r-robust, and demonstrated that the properties of r-minimum-degree, r- 
conne -ctivity and r-robustness share the same probability threshold in random /-partite graphs (see Fig. |l(b)[ ). 
Here, we will investigate a coarser rate of growth for the inter-layer edge formation p, and show that for such 
probability functions, the isoperimetric constant and the robustness parameter have the same asymptotic rate of 
growth. This will also play a role in the next section, where we investigate the algebraic connectivity of random 
interdependent networks. We start with the following lemma. 

Lemma 3: Consider a random /-partite graph G p = (V) U V9U • • • U V/,.• E p ) with node sets V) = { (i — l)n +1, (i — 
1 )n + 2,..., in] for 1 A i A k. Assume that the edge formation probability p satisfies limsup, woo ( k l ™ np A 1. 
Fix any e € (0, |]. There exists a constant a (that depends on p) such that the minimum degree d m .i n , maximum 
degree d max and isoperimetric constant i(G p ) a.a.s. satisfy 


oatvp A %(Gp) A d m i n A dmax A n(/ 1 )p 


^l + \/3 


/ In n \*~ e \ 
\(k-l)np) J 


( 6 ) 


The proof of the above lemma is given in Appendix |B] The above result leads to the following theorem. 
Theorem 2: Consider a random /-partite graph G with Bernoulli interconnections. Assume that the inter-layer 
edge formation probability p = pin) satisfies 


lim sup 

n—> 00 


In n 

(/ — 1 )np 


A 1. 


Then i(G) = 0(np) a.a.s., and furthermore, G is ©(np)-robust a.a.s. 

Proof: For inter-layer edge formation probabilities satisfying the given condition, we have i(G) = Q(np ) 
a.a.s. from Lemma [3] From Lemma [Q we have that the robustness parameter is bl(rip) a.a.s. Furthermore, since 
the robustness parameter is always less than the minimum degree of the graph, the robustness parameter is 0{np) 
a.a.s. from Lemma [3] which concludes the proof. ■ 

Once again, since adding edges to a network cannot decrease the isoperimetric constant or the robustness 
parameter, the above result immediately implies that for random interdependent networks with inter-layer edge 
formation probability satisfying limsup,,^^ ^-p)np < we ^ ave ? -( < ^) = - Z (np) and that the robustness is f l{np) 
a.a.s. This condition on the inter-layer edge formation probability has further implications for the structure of 
random interdependent networks. In next section, we use Lemma [3] to show that the algebraic connectivity of 
random interdependent networks scales as Q{np ) a.a.s. for all p that satisfy the condition, again regardless of the 
probability distributions over the topologies of the layers. 


V. Algebraic Connectivity of Random Interdependent Networks 

The algebraic connectivity of interdependent networks has started to receive attention in recent years. The 
authors of 0 analyzed the algebraic connectivity of deterministic interconnected networks with one-to-one weighted 
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symmetric inter-layer connections. The recent paper by Hernandez et al. studied the algebraic connectivity of a 
mean field model of interdependent networks where each layer has an identical structure, and the interconnections 
are all-to-all with appropriately chosen weights 0. Spectral properties of random interdependent networks (under 
the moniker of planted partition models ) have also been studied in research areas such as algorithms and machine 
learning fTTl . HU], H20H . ll23l . |j24]]. Here, we leverage our results from the previous section to provide a bound on 
the algebraic connectivity for random interdependent networks that, to the best of our knowledge, is the tightest 
known bound for the range of inter-layer edge formation probabilities that we consider. 

Theorem 3: Consider a random interdependent graph G = (G i, G'a,..., Gk,G p ) with Bernoulli interconnections 
and assume that the probability of inter-layer edge formation p satisfies limsup,^^ < 1- Then AafCi) = 

0 (np) a.a.s. 

Proof: In order to prove this theorem, we need to show that there exist constants 7 , (3 > 0 such that 7 np < 
X 2 (G) < (3np a.a.s. We start with proving the existence of constant 3. 

Consider the set of nodes V\ in the first layer. The number of edges between V\ and all other Vj, 2 < j < k is 
a binomial random variable B{n 2 {k — 1 ),p) and thus E[|<9TT|] = n 2 {k — 1 )p. By using the Chernoff bound ll25l 
for the random variable \dV\ |, we have (for 0 < 5 < 1) 

-E(|aVi |)« 2 

Pr(|<9TT| > (1 + <5)E[|c)TT|]) < e -3- . (7) 

Choosing 5 = —7=. the upper bound in the expression above becomes exp ^— w ■ Since Inn < n(k — 1 )p 

for n sufficiently large and for p satisfying the condition in the proposition, the right hand side of inequality ( 0 ) 
goes to zero as n —>• 00 . Thus \dV\ | < (1 + o(l))E[|c?Vi |] a.a.s. Therefore 

i(G) = min M < Mi < t 1 + °( 1 »" 2 (* ~ Dr, 

\A\<¥, \A\ | Ti| n 

ACV 1 L>V 2 U---UV k 

a.a.s. Using (0), we have the required result. 

Next, we prove the lower bound on A 2 (G). Consider the fc-partite subgraph of network G which is denoted by 
G p . By Lemma [3] and the inequality (0), we know that A 2 (G p ) > 7 np for some constant 7 asymptotically almost 
surely. Since adding edges to a graph does not decrease the algebraic connectivity of that graph ll26l . we have 
A 2 (G) > A 2 (G p ) > 7 np asymptotically almost surely. ■ 

Theorem 0 demonstrates the importance of inter-layer edges on the algebraic connectivity of the overall network 
when limsup^^oo (f}™ np < 1. This requirement on the growth rate of p cannot be reduced if one wishes to stay 
agnostic about the probability distributions over the topologies of the layers. Indeed, in the proof of the Theorem 
0 we showed that if lirri sup TWOO / k _pj np > 1 , a random fe-partite graph will have at least one isolated node 
asymptotically almost surely and thus has algebraic connectivity equal to zero asymptotically almost surely. In this 
case the quantity forms a coarse threshold for the algebraic connectivity being 0 , or growing as Q(np). 

On the other hand, if one had further information about the probability distributions over the layers, one could 
potentially relax the condition on p required in the above results. For instance, as mentioned in Section [TITJ when 
each of the ^-layers is an Erdos-Renyi graph formed with probability p, then the entire interdependent network is 
an Erdos-Renyi graph on kn nodes; in this case, the algebraic connectivity is il(np) asymptotically almost surely 
as long as lira sup n _ >0O < 1 |27t . This constraint on p differs by a factor of ^-y from the expression in 
Theorem 0 


VI. Summary and Future Work 

We studied the properties of r-robustness and algebraic connectivity in random interdependent networks. We 
started by analyzing random /.’-partite networks, and showed that r-robustness and r-minimum-degree (and r- 
connectivity) all share the same threshold function, despite the fact that r-robustness is a much stronger property 
than the others. This robustness carries over to random interdependent networks with arbitrary intra-layer topologies, 
and yields new insights into the structure of such networks (namely that they can tolerate the loss of a large number 
of nodes, and are resilient to misbehavior in certain dynamics). We also provided tight asymptotic rates of growth 
on the algebraic connectivity of random interdependent networks for certain ranges of inter-layer edge formation 
probabilities (again, regardless of the intra-layer topologies), showing the importance of the interdependencies 
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between networks in information diffusion dynamics. Our characterizations were built on a study of the isoperimetric 
constant of random interdependent and A;-partite graphs. 

There are various interesting avenues for future research, including a deeper investigation of the role of the 
intra-layer network topologies, and other probability distributions over the inter-layer edges (outside of Bernoulli 
interconnections). 


Appendix A 
Proof of Lemma [2] 

Proof: Suppose the edge formation probability is p = where x = o(lnlnn) and x -A oo 

when n -A oo. We have to show that for any set of vertices of size m, 1 < m < nk/2, there are at least m(r — 1) +1 
edges that leave the set a.a.s. 

Consider a set S C V\ U Vi U ■ ■ ■ U 14 with |5| = m. Assume that the set S contains s* nodes from V) for 
1 < i < k (i.e., IS 1 Cl V\ = Si > 0). Define Eg as the event that 777(7" — 1) or fewer edges leave S. Note that \dS\ 
is a binomial random variable with parameters f2i=i s i (St=i t^z( n — s t)l an< i P- 
We have that 


T: N ( n - s t> = 'Yh - 1) - m + Si) 

1=1 \t=l,t^l ) 1=1 


= n{k — 1)771 — m 2 + s h 


1=1 


( 8 ) 


Then we have 



n{k — 1)777 


m 


+ Eli * 


p 1 {i-p) 


— T d n ( fc - 1 ) m - m 2 +£w=i s i ! - i 


< 


< 



p\l — p) n ( A:_1 i m_m2+ ^'= lS *~ i 


p\l-p) 


/1 -1 \ (fc —l)m^ 

_ \n(fc—l)ra— -—77- 1 


(9) 


The first and the second inequality follow from the inequalities 0 < m 2 — Yli= 1 s z 2 < ^ P” 7 which is a 
straightforward result of the Cauchy-Schwartz inequality in ©. Next note that k > 2 and for I < i < m(r — 1) 
and sufficiently large n, we have 

(n(k-l-i 
^n(fe-l)m)y_i(i _ p'y l {k—l)m— V*- 1 )™ 2 —(i—1) 

n(k — l)m — i + 1 p 

= — --- X - 

i 1 — p 

n{k — l)m — m(r — 1) + 1 p 

— - 7 -r-x- 

777(7" — 1 ) 1 — p 

n(k — 1) — r + 1 p 

> — ---x- 

7" — 1 1 — p 

> 1 . 


Thus there must exist some constant R > 0 such that 

m(r— 1 ) 


11 J. J * , , >. 

Pr(E s )< £ 

1=0 ^ * ' 
n(k - l)m\ (r _i) 
m(r — 1) ) P 1 11 


\n(k—l)m— —— 


( 10 ) 


—m(r— 1 ) 
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Define P m as the probability that there exists a set of nodes T such that \T\ = m and \ 0 T\ < m(r — 1 ). Then 
using the inequality (™J < (^-) m yields 

Pm < J2 Pr (^) 


|5|=m, 

S CXJ k i=1 Vi 


< R 

< R 

= R 

< R 


nk\ (n(k — 1 )m 
mj \ m(r — 1) 


P 


m(r— 1) 


(i -p) 


n(k—l)m ———- in{r— 1) 


{nke\ m ( n{k - l)mep \ m{r 1} _ p )n(fc-i)m-^=^-m(r-i) 

\ m ) \ m(r — 1) / 

ke r n(l — p) n(fe_1) (n{k — l)p) r_1 


( 11 ) 


r—1 \ m 


(r - l) r - 1 (l -py - 1 

c\n(l — p) n ( fc_1 ) (n(k — 1 )p) rl 
m{ 1 — p) ( k 


,, sPLzl > 

m(l — p) <« 


In the last step of the above inequality, ci is some constant satisfying 

/ce r 2fce r 


(r — l) r_1 (l — p) 


—i < ci < 


(r- 1) 


r—1 ’ 


_ lnn+(r—1) lnln n+x 


for sufficiently large n. Recalling the function pin) = - (k-i)n 

yields 

/ cine~ n( ' k ~ 1 ' lp (n(k — l)p) r_lN m 

< xt -- 


and using the inequality 1 — p < e 


-p 




= R Ci 


(fc-l)m 

m(l — p) k 
In n + (r — 1) In In n + x 


In n 

\ m 


r—1 


(fc-l)m 

m(l — p) k 


< R 


c 2 e“ 


my 1 — p) fc 


Due to the fact that 

Inn, 
r—1 


< 2 for sufficiently large n, we can consider a constant upper bound 

C2 for ci ( In n+ ^ r ^^ n ln In n+X \ such that 0 < c 2 < ci 2 r_1 . Next, we substitute the Taylor series expansion 

ln(l — p) = — yyy I y for p G [0, l) in the above inequality to obtain 


Pm. < R 


= R 




m 

(fc — 1 )mp 


c 2 e x e >= 


e*p{^r 2 £“ 2 GT 


771. 


Since we have £“ 2 C < £“ 2 p- 2 = + and ^G^p 2 


0 as n —>• oo, there must exist a constant C3 such 


that 0 < 


< -—pP-p 2 Y^S=2 < C3 < 1 for sufficiently large n. Therefore, 


Pm < R \ C 2 e' 


C3 


_ (fc-l)mp \ m 

e e k 


m 


= R C4- 


(fc-l)mp \ m 

e~ x e fc \ 


771. 


where 0 < C4 = c 2 e C3 < 


fce’’ +1 2 r ' 


(i 1) T ’— 1 ' 

Consider the function f{m ) = 


(fe-l)mp 

- Since = 

m dm 


df _ 


and ^ > 0 for 777 > Therefore, f(m ) < max{/(l), /(nfc/2)} for m G {1,2,..., Ln,A:/2j }. We have 


, we have that < 0 for m < 


(fc —l)nfcp 

f(nk/ 2 ) = — 2I ~ = ^g(- ln «+(r’- 1 ) lnln «+a:)/ 2 _ 

nk k 
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/. \ „ /x__ (fc-l)p 

Since In Inn = o(lnn), we have that f(nk/ 2) = o(l). Moreover, 1 < /(1) = e * < e and thus for sufficiently 

large n we must have f(m) < /(1) < e. Therefore, 

Pm < R{c4e 1 ~ X ) m . 


Let P 0 be the probability that there exists a set S with size nk/2 or less that it is not ?’-reachable. Then by the 
union bound we have 


\nk/2\ oo 

p 0 < J2 ^ E p( c ^~ x r 

m =1 m= 1 


Rc A e'- x 
1 - c 4 e 1_3: 


o(l), 


since i -too as n -> oo. Thus we have the required result. 


Appendix B 
Proof of Lemma [3] 

Proof: The inequality i(G p ) < d m i n is clear from the definition of the i(G p ). We will show that d max < 
n{k — 1 )p ^1 + \/3( (fcl[y) n p )^~ £ ) asymptotically almost surely. Let dj denote the degree of vertex j, 1 < j < kn. 
From the definition, dj is a binomial random variable with parameters n(k — 1) and p and thus E [dj] = n(k — 1 )p. 
Then, for any 0 < (3 < y/3 we have li25ll . ll27l 

E[d,]/3 2 

Pr (dj > (1 + P)E [dj]) < e -3—. (12) 

Choose P = \/?>( (i^\y np ) 1 -' i which is at most y/3 for p satisfying the conditions in the lemma and for sufficiently 
large n. Substituting into equation (fl2l) . we have 

Pr (dj > (1 + p)E[dj]) < e " (lnn)( W=WT )_2e . 

The probability that d max is higher than (1 + /3)E [dj] equals the probability that at least one of the vertices has 
degree higher than (1 + /?)E [dj], which by the union bound is upper bounded by 


Pr (dmax > (1 + P)E[dj]) < knPr(dj > (1 +/3)E [dj]) 

< fc e (lnn)_(lnn)( A^TT ) " 26 


Since the right hand-side of the above inequality goes to zero as n —» oo for p satisfying the condition in the 
lemma, we conclude that 

Inn 


dn 


< n(fc - l)p 1 + >/3 


(k — 1 )np 


asymptotically almost surely. 

Next we aim to prove the lower-bound for i(G p ) in ©. In order to prove this, we show that for any set of 
vertices of size m, 1 < m < nk/2, there are at least amnp edges that leave the set, for some constant a that we 
will specify later and probability p satisfying linq^oc {k-V)np < ■*“ 

Consider a set S' C Vi U V 2 U ■ ■ ■ U 14 with |/S'! = m. Assume that the set S contains nodes from V) for 
1 < i < k (i.e., \S Cl Vi\ = s* > 0). Define Eg as the event that amnp or fewer edges leave S. Note that \dS\ is 

a binomial random variable with parameters Y^i=i s i (lEt=\ t^i( n — j and p. Similarly to the equality © and 
inequality @, we have that 



P*(l ~P) 




(fc-1)" 


(13) 
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Next note that k > 2 and for 1 < i < [amnp ], 

(n(k-l_ p^n(k—l)m— (fc ~* )m2 -i 


_ p^n(k-l)m- (k * )m2 -(a-1) 

n(k — l)m — i + 1 » 

= — : - x 


1 -P 


n(& — l)m — amnp +1 p 

^-x 

amnp 

k — 1 — an 1 

> -- x 


1 — p 


a 


1 — p 


1 

> 

a 


for a < 1 which will be satisfied by our eventual choice for a. 

Now let P m denote the probability of the event that there exists a set of size m with [amnp\ or fewer number 
of edges leaving it. Then there must exist some constant R > 0 such that by the same procedure as in inequalities 
Cl and CD. we have 


Pm < R 


< R 


nke\ m ( n(k — l)mep \ amnp 
m J \ amnp ) 

e m ln ( ) e ~p(n(k-l)m-amnp- ) 


(k- l) e ^ amnp 


a 


where 


= Re mh ^ m \ 

h{m) = anp + anp\n{k — 1) — anplna + In 

/ 


(14) 


nke\ (k — 1 )m, 

— p[n[k — 1) — anp — 


m 


k 


, , (k — 1 )pm . 

= 1 + In k H-In m + np 

k 


\ 


Inn 


a + a ln(fc — 1 ) — a In a H- (k — 1 ) + ap 

np 

-v- 

r(a) 


V 


(15) 




Since 


dhPpJ = (k i) P _ X i s negative for m < ( k -i)p and P ositive for m > (k-i)p ’ we have 
h(m) < max{fi(l), h(nk/ 2 )} 


< max 11 + In k + ^ + npT(a), 1 + In 2 + np ^T(a) + 


(k — 1) Inn 


in n \ 
np J ) 


From (flSb . = ln(A; — 1) — In a +p > 0 and thus T(a) is an increasing function in a for a < (k — 1), with 

T(0) = ^ — (k — 1) which is negative and bounded away from 0 for sufficiently large n (by the assumption on 
p in the lemma). Therefore, there exists some sufficiently small positive constant a such that h(rri) < —anp for 
sufficiently large n. Thus CD becomes P m < R e ~a mn P for sufficiently large n. 

The probability that i(G p ) < anp is upper bounded by the sum of the probabilities P rn for 1 < m < nk /2. 
Using the above inequality, we have 

nk /2 nk/2 


Pr (i(G p ) < anp) < 


m= 1 


m= 1 
oo 


< 


P E 


3 —amnp 


amnp 


m= 1 


= R- 


0 —anp 


^ _ g—cunp 7 


which goes to 0 as n —» oo. Therefore, we have i(G p ) > anp asymptotically almost surely. 
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